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Abstract 

The non-Gaussian measures of the particle number fluctuations - skewness Sa and kurtosis na 2 
- are calculated in a vicinity of the critical point. This point corresponds to the end point of the 
first-order liquid-gas phase transition. The gaseous phase is characterized by the positive values of 
skewness while the liquid phase has negative skew. The kurtosis appears to be significantly negative at 
the critical density and supercritical temperatures. The skewness and kurtosis diverge at the critical 
point. The classical van der Waals equation of state in the grand canonical ensemble formulation is 
used in our studies. Neglecting effects of the quantum statistics we succeed to obtain the analytical 
expressions for the rich structures of the skewness and kurtosis in a wide region around the critical 
point. These results have universal form, i.e., they do not depend on particular values of the van der 
Waals parameters a and b. The strongly intensive measures of particle number and energy fluctuations 
are also considered and show singular behavior in the vicinity of the critical point. 
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I. INTRODUCTION 


The first-order phase transition is among the most general phenomena in physics. This 
phase transition exists in atomic and/or molecular systems, in the system of interacting nu¬ 
cleons (nuclear matter), and, most probably, in the QCD, between hadrons and quark-gluon 
plasma at large baryonic densities. The liquid-gas phase transition line in the plane of tem¬ 
perature T and chemical potential /i has the end point, which is called the critical point (CP). 
The CP demonstrates some universal features typical for the second-order phase transitions, 
particularly, anomalously large fluctuations. The study of event-by-event fluctuations in high- 
energy nucleus-nucleus collisions opens new possibilities to investigate properties of strongly 
interacting matter (see, e.g., Refs, jl] and [ 2 ] and references therein) and the experimental 


search for the QCD CP is now in progress (see, e.g., Re 
fluctuation signals of the QCD CP were discussed in Ref. 


a, 


and references therein). The 


and higher moments of fluctua¬ 
tions of conserved charges were suggested as probes to study the phase structure of QCD { 5 , 
Particularly, the higher-order non-Gaussian measures such as the skewness Sa and kurtosis 


na 2 of conserved charges fluctuations have attracted much attention (see, e.g., Ref 


'. [7] and jg 


8 ]). 


Experimentally, the STAR collaboration has measured the higher moments of net-proton and 
net-charge multiplicity in Au+Au collisions {9 11] . See also recent review 12] on the search 
for critical behavior of strongly interacting matter at the CERN Super Proton Synchrotron. 


Calculations of 
QCD models 


righer moments of conserved charges has been performed in various effective 


13H17I] . The effects of non-equilibrium evolution of these observables in heavy-ion 


collisions have also been considered recently 18]. Another fluctuation measures of interest are 
the so-called strongly intensive quantities jl9[ which are normally not sensitive to the fluctu¬ 
ations of the system volume. This is especially relevant for heavy-ion collision experiments 
where size of the colliding system varies strongly on event-by-event level. 

The van der Waals (VDW) equation is a simple analytical model of the pressure function 
p for equilibrium systems of particles with both attractive and repulsive interactions. In the 


canonical ensemble (CE) it reads as (see, e.g., Refs. 20j, 2l|), 


p(T,n) = 


NT 


N 2 


n T 


V-bN a V 2 


— an 


( 1 ) 


1 — bn 

where n = N/V is the particle number density while the VDW parameters a > 0 and b > 0 
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describe the attractive and repulsive interactions, respectively. The first term on the right- 
hand-side of Eq. (PP) corresponds to the excluded volume (EV) correction, which manifests 
itself in a substitution of a total volume V by the available volume, V av = V — b N. The second 
term comes from the mean held which describes attractive interactions between particles. With 
regards to its asymptotic behavior in the vicinity of the critical point the VDW model belongs 
to the mean-held theory universality class. The number of particles N is hxed in the CE. In 
order to apply the VDW equation of state to systems with variable number of particles and 


calculate their huctuations the grand canonical ensemble (GCE) formulation is needed, 
procedure was hrstly performed for the EV model, i.e., for a = 0 in Eq. ([]]), in Refs. [22 


’his 


23]. 


In our recent paper 24j> the full VDW equation (JT[) , with both attractive and repulsive terms, 
was transformed from the CE to the GCE for systems with Boltzmann statistics while the 
formulation which properly includes effects of quantum statistics was obtained in Ref. 25]. 


Note that the EV and VDW models can also be conveniently treated within the GCE in a 


framework of the thermodynamic mean-held approach (see Refs. 26-28])^ 


In the present paper we use our recent results of the GCE formulation 24] as a starting point 
for calculating the scaled variance, skewness, and kurtosis of particle number huctuations, and 
also the strongly intensive measures of particle number and energy huctuations. The Boltzmann 
approximation is adopted. This gives a possibility to obtain the analytical expressions for the 
universal structure of huctuations in a vicinity of the CP. 

The paper is organized as follows. In Sec. [IT] the VDW equation of state is formulated within 
the GCE. In Sec. IIIII the particle number huctuations - scaled variance, skewness, and kurtosis - 
are calculated and their behavior in a vicinity of the CP is analyzed. In Sec. IIVI the strongly 
intensive measures of huctuations for the energy and number of particles are considered. A 
summary in Sec. |V] closes the article. 


II. VDW EQUATION IN THE GCE 

The canonical ensemble VDW pressure function in (1) corresponds to the Boltzmann ap¬ 
proximation. This expression is used in the present paper, i.e., effects of the quantum statistics 
(Bose or Fermi) will be neglected. The VDW pressure is a unique function of variables T and 
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n for all T > 0 and 0 < n < 1/b, and this equation of state contains the first-order liquid-gas 
phase transition and has the CP. The CP in (T, n)-plane, i.e. the point (T c ,n c ), corresponds to 
the temperature and particle number density, where the following derivatives are equal to zero, 

d 2 p 


A) = 0 

dn ) T 


dn 2 


= 0 . 


and the thermodynamical quantities at the CP are equal to: 


T r = 


8 a 


n, = 


Pc = 


27 b 2 


27 b 1 L 3 b 

At T > T c the following inequality is always valid, 

g) >o. 

an ) T 

while at T <T C there appears an unstable interval [ 71 . 1 , 77 . 2 ] with 

dp' 


dn 


< 0 . 


( 2 ) 


(3) 


(4) 


(5) 


This means that the VDW isotherm p(T, n) has a local maximum at n = 77.1 and a local 
minimum at n = 77.2 > 72 X for T < T c . The unstable part (J5J) of the VDW isotherm at the interval 
[77.1,772], together with two additional metastable parts - [?7 s ,77i] and [t7 2 ,77/] - are transformed 
to a mixture of two phases: a gas with density n g < n\ and a liquid with density ni > 772 . This 
is done in accordance with the Maxwell rule of the equal areas (see, e.g., Refs. [20, |21]) which 
leads to a constant pressure p(T,n g ) = p(T,n{) inside the density interval [t7 9 ,77;]. 

In the GCE the pressure should be defined in terms of its natural variables: temperature 
T and chemical potential /j. The p(T, p) function contains a complete information about 
equilibrium physical systems. Other thermodynamical quantities, such as particle number 
density n(T,p), entropy density s(T,p), and energy density e(T, / 1 ) can be presented in terms 
of p and its T- and yU-derivatives: 








dp\ 

+ li{ Mdr~ P - 


( 6 ) 


The VDW equation of state in the GCE is obtained in the form of transcendental equation 
for particle number density n = n(T,p ) as a function of T and p (see Ref. 24] for details): 

n ld (T,p*) „ „ bn 


n(T,p ) = 


1 + bn id (T,p*) 


p* = p - T 


1 — bn 


+ 2(777 , 


(7) 
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where n ld is a particle number density in the ideal Boltzmann gas 


n id (T,/i) = exp(^ 


p \ dm 2 T 


T) 2vr 2 




m 


( 8 ) 


with d being the degeneracy factor and m being the particle mass. The K 2 (x) is the modi¬ 
fied Bessel function of the second kind. The GCE VDW pressure p(T, /i) is then obtained by 
inserting n(T,p) from (J7J) into Eq. (JT]). Note, the relativistic form of a dispersion relation is 
considered, ui(k) = \Jm 2 + k 2 , where u and k are the free single-particle energy and momen¬ 
tum, respectively 1 . This makes the present formulation being suitable for high-energy physics 
applications. 

In the GCE there is a unique solution of Eq. ((7]) at T > T c , while at T < T c it may have 
either one solution or three different solutions for particle number density n(T, p). In that 
case the solution which corresponds to a largest pressure survives in accordance to the Gibbs 
criterion. The liquid-gas mixed phase in the T-p plane belongs to the line p = p c {T), where 
two solutions with different particle number densities, n g (T,p) and ni(T,p), correspond to the 
equal pressures, p g (T, p) = pi(T, p). The Maxwell rule of the equal areas and the Gibbs criteria 
of equal pressures for the gas and liquid at the same T and p values appear to be the equivalent 
descriptions of the first-order liquid-gas phase transitions (see Ref. [25] for details). 


III. PARTICLE NUMBER FLUCTUATIONS 

A. Fluctuations in the GCE 

Let the particle number A be a random variable with the normalized probability distribution 
V(N). The fc-th moment ( N k ) is then defined as 

(N k ) = ^N k V(N) . (9) 

N 

Let us introduce the variance, a 2 = ((AiV) 2 ), where AN = N — ( N ). The scaled variance, 

" (V) ■ (10) 


1 We use the system of units, where the Plank constant H , the speed of light c, and the Boltzmann constant k 
are equal to unity, h = c = k = 1. 
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characterizes the width of the V(N) distribution. Note that oj[N] = 1 for the Poisson distribu¬ 
tion V(N) = exp(-(iV)) (N) n /N\. 

The skewness Sa is defined as 


S(T = 


((A AO 3 ) 


a* 


( 11 ) 


The skewness measures the degree of asymmetry of the distribution V(N) around its mean 
value (N). Positive skewness indicates a distribution with an asymmetric tail extending more 
to the right , i.e., toward iV-values with N > (N). Negative skewness indicates a distribution 
with an asymmetric tail extending more to the left, i.e., toward iV-values with N < (N). If the 
V(N) distribution is symmetric around its mean value, i.e., the right and left tails are equal, it 
has zero skewness. This is the case for the normal Gaussian distribution, whereas the Poisson 
distribution shows a positive value of the skewness, Sa = 1. 

The (excess) kurtosis na 2 is the measure of “peakedness” of the probability distribution 
V(N), 


2\2 


na = 


((AAO 4 ) - 3 ((AN) 2 ) 


a A 


( 12 ) 


The kurtosis (H3 measures the degree to which a distribution is more or less peaked than 
a normal Gaussian distribution. Positive kurtosis indicates a relatively peaked distribution. 
Negative kurtosis indicates a relatively flat distribution. For the Poisson distribution one has 
positive value of the kurtosis, na 2 = 1. 

The normal Gaussian distribution corresponds to the zero value of both the skewness (fTTl) 
and the (excess) kurtosis na- Therefore, (strong) deviations of Sa and/or na 2 from zero are 
the signatures of the (highly) non-Gaussian shape of the particle number distribution V(N). 

In the GCE the system is defined by the pressure p given in terms of its natural variables T 
and pi. The particle number fluctuations can be characterized by the following dimensionless 
cumnlants (susceptibilities), 

d n (p/T i ) 


kn d(p/T) n ’ 

which are connected to the moments of the particle number distribution as 


h = 


(N) 
VT 3 ’ 


ko = 


((AAQ 2 ) 

VT 3 


k?. = 


((AAQ 3 ) 

VT 3 


2\2 


hi = 


((AAO 4 ) - 3((AA0 2 ) 

VT 3 


(13) 


(14) 
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where (...) denotes the GCE averaging. The scaled variance (fTOlh skewness cm. and kurtosis 
(fl2l) are the intensive fluctuation measures that remain finite in the thermodynamic limit V —> 
oo. They can be expressed in terms of the susceptibilities as the following 


cu[N] = 


k 2 
h ’ 


S (j = 


h 

k 2 


na 2 = 


(15) 


B. Scaled variance 

The scaled variance defined in Eq. (1T5]) can be calculated through the /i-derivative of the 
particle density. For pure phases in the classical VDW gas it results in (see also Ref. 24]) 

T , ^ 

(16) 


(17) 


/ dn\ 

1 

2 an 

\dp) T 

(1 — bn) 2 

T 


n 

In terms of the reduced quantities, T = T/T c and n = n/n c , Eq. flTfil) reads as 

-l 


w [iV] = - 


n 

AT 


9 [(3 — n) 2 

and possesses the universal form independent of the specific values of the VDW parameters a 
and b. From Eq. (1171) it follows that u;[iV] —> 1 at n —> 0 (this corresponds to the ideal gas 
limit and the Poisson V(N) distribution), and cu[iV] —>■ 0 at n —)■ 3 (this corresponds to the 
liquid with the highest possible density). According to its definition, the scaled variance u;[iV] 
is a positive quantity. This is indeed the case for all T and n values that correspond to stable 
and even metastable states. The scaled variance diverges, oj[N] oo, at the CP. Introducing 
quantities p — n — 1 and r = T — 1 one finds at r C 1 and p 1: 


= 9 


t+ 4 P +Tf> 


-1 


(18) 


The scaled variance (TT6l) as a function of T and n , is plotted in Fig. [T] for both stable and 
metastable pure phases. 

C. Skewness 


The skewness Sa can be calculated as 


= HN}f 


1 — 3 bn 
(1 — bn) 3 


(19) 
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Figure 1: (Color online) The scaled variance u;[iV] (1161) on the ( n,T ) phase diagram for both stable 
and metastable pure phases. Several lines of constant values of cu[lV] are shown. The grey area depicts 
region where pure phase is mechanically unstable. 


In the reduced variables it reads 


S<r = - 


1 n 

-2 

1 — n 

(3 — n) 2 4 f_ 


(3 — n ) 3 


( 20 ) 


As a function of the reduced temperature and density the skewness Sa is plotted in Fig. [ 2 ] 
for both stable and metastable pure phases. It is clearly seen from Eqs. (TIU1) and (12Uj) that 
the skewness is positive at n < 1 (the gaseous phase), negative at n > 1 (the liquid phase), 
and Sa — 0 at n — 1. At n —> 0 one finds that Sa —>■ 1 . This is a small asymmetry of the 
particle number distribution and it corresponds to the Poisson distribution that takes place for 
the ideal Boltzmann gas. 

In a vicinity of the CP one finds, 

-2 


Sa ^ 


;P 


T + ^P + T P 


( 21 ) 


The CP, T — n — 1, i.e., p = t = 0, is a point of the essential singularity of the skewness 
measure. For example, at r = 0 one finds from Eq. (1211) that Sa —l + oo at /) ->• + 0 and 
Sa —> — oo at p —> — 0. At the same time, Sa = 0 at any r > 0 and p = 0. 
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Figure 2: (Color online) The same as in Fig. Q] but for the skewness Sa (1201) . Different values of 
skewness are illustrated by the typical V(N) distributions on the right panel. 

D. Kurtosis 

The kurtosis ncr 2 can be calculated as 

“ 2 = | = ^> 2 + T (fr) r = (22) 

As a function of the reduced temperature and density the kurtosis na 2 is plotted in Fig. [3] for 
both stable and metastable pure phases. It is seen from these figures that at T < 1 the kurtosis 
is positive (leptokurtic) for both n < 1 (the gaseous phase) and n > 1 (the liquid phases). 
Approaching the CP the kurtosis diverges. At r = 0 and p ■ Cl one finds 

k < j 2 oc p -6 . (23) 

Notably, the kurtosis has attains large negative values (platykurtic) at critical density n — 1 
and temperatures just above the critical, T > 1. In this region, one finds 

ko 1 oc - r -3 . (24) 

This indicates that particle number distribution has a flat peak in that region, much flatter than 
the corresponding Gaussian with the same width. This region can be identified as the crossover 
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Figure 3: (Color online) The same as in Figs. [T] and [2] but for the kurtosis k<j 2 (1221) . Different values 
of kurtosis are illustrated by the typical V(N) distributions on the right panel. 


region, where rapid although smooth transition between gaseous and liquid phases takes place. 
This is in line with arguments that the crossover region near the CP is characterized by the 
negative sign of kurtosis Q. In the vicinity of the CP the kurtosis changes rapidly and attains 
both positive and negative values. 

At n —» 0 the VDW equation of state corresponds to the ideal Boltzmann gas limit. In this 
case, V(N) approaches the Poisson distribution, and, as it follows from Eqs. (1201) and (122] ) . 
the skewness and kurtosis both approach the Poisson expectation values, i.e., Sa —> 1 and 
na 2 —y 1. It should be noted that the ideal Boltzmann gas with the Poisson distribution reveals 
small deviations from the Gaussian V(N) distribution for which, by construction, Sa = 0 and 
na 2 = 0. 


The nucleon number fluctuations were recently considered in Ref. 29] within the modified 


VDW equation. To describe the nuclear matter ground state the Fermi statistics of nucleons 
was introduced 251 ]. The Boltzmann approximation is assumed in the present study and it 


corresponds to the original version of the VDW equation. This gives a possibility to obtain 
the analytical expressions for the scaled variance, skewness, and kurtosis. Note also that our 
results for the particle number fluctuations presented in Figs. |T]{3] are universal, i.e., they are 
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independent on the specific numerical values of the VDW parameters a and b. Therefore, they 
can be applied to very different physical systems - from the CP of water with T c = 647 K° up 
to the CP in of the nuclear matter with T c = 2 x 10 12 K°. The universality of the results is lost 
if one would take into account the effects of quantum statistics. On the other hand, in most 
cases the inclusion of Fermi statistics changes the results on a quantitative level, qualitatively 
they remain essentially the same as in case of the Boltzmann statistics. 


IV. STRONGLY INTENSIVE QUANTITIES 


The results of the previous section demonstrate a strong increase of the particle number 
fluctuations in a vicinity of the CP. The fluctuations may become also very large for metastable 
states, i.e., super-cooled gas and/or super-heated liquid. These fluctuation signals from phase 
transitions in the nuclear matter can be observed in the event-by-event analysis of heavy ion 
collisions. Note, however, that in these processes there is one more source of particle number 
fluctuations, namely, event-by-event fluctuations of the system volume. These volume fluctua¬ 
tions are mainly of the geometrical origin, and they can hardly be avoided in nucleus-nucleus 
reactions. Thus, one may observe large experimental fluctuations even for simple statistical 
systems, e.g., the ideal gas. The strongly intensive measures of the fluctuations defined in 


terms of two extensive quantities A and B were suggested in Ref. 


19]. For statistical systems 


in a case of the absence of phase transitions these measures within the GCE formulation are 
independent of the system volume and its fluctuations. Note, however, that in systems with 
the CP (in general, for the 2nd order phase transitions) the critical behavior does depend on 
the system volume and shows the characteristic finite-size scaling. This implies that strongly 
intensive quantities are also volume-dependent near the CP. Thus, using the strongly intensive 
measures one excludes trivial volume fluctuations for normal statistical systems, and a presence 
of large fluctuations in terms of these measures can be considered as an indication of critical 
behavior. 

In the present paper we consider the strongly intensive measures of total energy E and 
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particle number N fluctuations for the VDW equation of state. They are defined as 


A[E,N] = Q 1 (JV)u[E] - (E)u[JV] 

E[B,JV] = Cp[{iV)w(£;] + (£)w[lV] - 2[(EN) - {E) (N) 


(25) 

(26) 


where C A and C^ 1 are the normalization factors that have been suggested in the following 


form [30] 


C A = C E = (iV)w[e] . 


(27) 


with u[e\ being the scaled variance of a single-particle energy distribution in the VDW system. 
To proceed it is necessary to calculate uj[e], uj[E], and ( EN). 

In the VDW gas the average single-particle energy e is independent of the parameter b, but 
it is modified due to a presence of the attractive mean field: 

- - frr\ N o rr , K l (m/T) 

e = e id T -a— = 3 T + m——— - an, 28 

V A 2 (m/T) 

where Td is the average single-particle energy in the relativistic ideal gas. The variance of the 
single-particle energy is insensitive to the presence of the VDW mean field, and one obtains 


w[e] 



I dT 


(29) 


The mean total energy is 

(E) = ( (V - a*p) N^ = e u (N) - £ (N 2 ) = i id (N) - £ (/V) 2 - a ^ ~ < JV >' . (30) 

The first and second terms in the right hand side of Eq. (I30p are proportional to ( N ). On the 
other hand, the third term remains finite outside the critical point in the thermodynamic limit 
V — y oo. Therefore, one obtains 


(E) = (e id - an) ( N ) . 


For uj[E\ one then finds 

(A 2 ) - (A) 2 


W [A] = 


(E) 


J_ T 2 (d{E)\ 

< E ) V dT )^ /T 


w[e] + 


(e id - ‘Ian) 2 
6 id — an 


w[iV] . 


(31) 


(32) 
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Finally, the correlations between E and N can be calculated as the following 


(• EN) - (E) (N) = T 2 = (e.d - 2on) ( N) w[JV] . (33) 

Substituting the above formulae into Eqs. (125]) and (126]) one finds the following expressions 
for the strongly intensive quantities: 


A[ g ,iV] = 1 - (34) 

4 - 4 

E[E,fV] = 1 + == — uj[N] . (35) 

^id e id 

In the absence of the attractive interactions (i.e., a = 0), one can readily see from Eqs. (j34|) and 
(1351) that A [E,N] = £[E,iV] = 1, thus, in the excluded-volume model the strongly intensive 
quantities are the same as in the ideal Boltzmann gas. 




n 

Figure 4: (Color online) The strongly intensive quantity A [E, N] (1361) on the (h, T) phase diagram for 
both stable and metastable pure phases. Several lines of constant values of A [E,N] are shown. 

The expressions (l34j) and (l35|) for A[E, A] and T,[E,N] become more transparent if one 
considers the non-relativistic limit, 6;d = 3T/2 and ef d — ef d = 3T 2 /2. Note that the rest energy, 
m, has been subtracted, thus, only the kinetic energy fluctuations contribute to e 2 d . One then 
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Figure 5: (Color online) The strongly intensive quantity T,[E,N] (1371) on the ( n,T ) phase diagram for 
both stable and metastable pure phases. Several lines of constant values of T,[E,N] are shown. 


obtains 


A [E,N] = 1 


2 an(3T — 3 an) 

3 


cu[N] = 1 


9 n 
AT 


9 n 

8 f 




2 (1^71^ 27 77^ 

E[£,iV] = l + -— w [iV] = l + --=-u{N] 


(36) 

(37) 


The quantities A [A, A] and E[E,N] are depicted in Figs. [I] and EJ respectively. Both the 
T[E,N] and A[E,N] measures approach unity in both zero density, n —> 0, and packing, 
n —)■ 3, limits and diverge at the CP. Note that the measure is always positive and 

even larger than unity, while the A[E, A] measure attains both positive and negative values. 


V. SUMMARY 

Particle number fluctuations up to the fourth order - scaled variance cu[iV], skewness S'cr, and 
kurtosis na 2 - have been calculated for the classical van der Waals equation of state. Analytical 
formulae were obtained and used for an analysis of the fluctuation behavior in a vicinity of the 
CP. The obtained results have the universal form, i.e., they are independent of the specific 
numerical values of the VDW parameters a and b. 
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The skewness Sa is positive at n < n c and negative at n > n c for all values of temperature. 
This means that gaseous and liquid phases at T <T C are clearly characterized by positive and 
negative values of Sa, respectively. Above the critical temperature T c the skewness is equal 
to zero at the n = n c line, and this line can be associated with the crossover transition from 
gaseous to liquid matter. 

The kurtosis ko 1 2 3 4 is very sensitive to the proximity of the critical point and also has a rich 
structure. At T < T c the kurtosis is positive in both phases. In the vicinity of the CP the 
kurtosis changes rapidly and can attain both positive and negative values at T > T c . Notably, 
in the crossover region just above T c it attains large negative values. 

The strongly intensive measures A[E,N] and T,[E,N] for the fluctuations of system energy 
and number of particles have been also calculated. In the excluded-volume model, i.e. at a = 0, 
these measures are the same as in the ideal gas and are equal to unity. In the full van der Waals 
equation, however, they diverge at the critical point and differ significantly from unity in the 
region around it. These measures can be used to study the fluctuations in the nuclear matter 
created in heavy ion collisions. 

Acknowledgments 

We are thankful to M. Gazdzicki for fruitful comments and discussions. This work was 
supported by the Humboldt Foundation, by the Program of Fundamental Research of the 
Department of Physics and Astronomy of National Academy of Sciences of Ukraine, and by 
HIC for FAIR within the LOEWE program of the State of Hesse. 


[1] V. Koch, in Relativistic Heavy Ion Physics, Landold-Bornstein Volume 1/23, edited by R. Stock 
(Springer, Berlin, 2010). 

[2] M. I. Gorenstein, arXiv: 1505.04135 [nucl-th]; 9th International Workshop CPOD2014. 

[3] M. Gazdzicki, M. I. Gorenstein, and P. Seyboth, Int. Journ. Mod. Phys. E 23, 1430008 (2014). 

[4] M. Stephanov, K. Rajagopal, and E. Shuryak, Phys. Rev. Lett. 81, 4816 (1998); Phys. Rev. D 
60, 114028 (1999); M. Stephanov, Acta Phys. Polon. B 35, 2939 (2004). 


15 



[5] S. Ejiri, F. Karsch, and K. Redlich, Phys. Lett. B 633, 275 (2006). 

[6] M. Asakawa, S. Ejiri, and M. Kitazawa, Phys. Rev. Lett. 103 , 262301 (2009). 

[7] M. A. Stephanov, Phys. Rev. Lett. 102 , 032301 (2009). 

[8] M. A. Stephanov, Phys. Rev. Lett. 107, 052301 (2011). 

[9] M. M. Aggarwal et al. [STAR Collaboration], Phys. Rev. Lett. 105 , 022302 (2010). 

[10] L. Adamczyk et al. [STAR Collaboration], Phys. Rev. Lett. 112 , 032302 (2014). 

[11] L. Adamczyk et al. [STAR Collaboration], Phys. Rev. Lett. 113 , 092301 (2014). 

[12] M. Gazdzicki and P. Seyboth, arXiv: 1506.08141 [nucl-ex], 

[13] B. J. Schaefer and M. Wagner, Phys. Rev. D 85, 034027 (2012). 

[14] K. Morita, B. Friman, K. Redlich, and V. Skokov, Phys. Rev. C 88, 034903 (2013). 

[15] K. Morita, V. Skokov, B. Friman, and K. Redlich, Eur. Phys. J. C 74, 2706 (2014). 

[16] K. Morita, B. Friman, and K. Redlich, Phys. Lett. B 741, 178 (2015). 

[17] S. Chatterjee and K. A. Mohan, arXiv: 1502.00648 [nucl-th], 

[18] S. Mukherjee, R. Venugopalan, and Y. Yin, Phys. Rev. C 92 , 034912 (2015). 

[19] M. I. Gorenstein and M. Gazdzicki, Phys. Rev. C 84, 014904 (2011). 

[20] W. Greiner, L. Neise, and H. Stocker, Thermodynamics and Statistical Mechanics , 1995 Springer- 
Verlag New York, Inc. 

[21] L. D. Landau and E. M. Lifshitz, Statistical Physics (Oxford: Pergamon) 1975. 

[22] M. I. Gorenstein, V. K. Petrov, and G. M. Zinovjev, Phys. Lett. B 106 , 327 (1981). 

[23] D. H. Rischke, M. I. Gorenstein, H. Stoecker, and W. Greiner, Z. Phys. C 51 , 485 (1991). 

[24] V. Vovchenko, D. V. Anchishkin, and M. I. Gorenstein, J. Phys. A 48, 305001 (2015). 

[25] V. Vovchenko, D. V. Anchishkin, and M. I. Gorenstein, Phys. Rev. C 91 , 064314 (2015). 

[26] D.V. Anchishkin, Sov. Phys. JETP 75, 195 (1992) [Zh. Eksp. Teor. Fiz. 102, 369 (1992)]. 

[27] D. Anchishkin, E. Suhonen, Nucl. Phys. A 586, 734 (1995). 

[28] D. Anchishkin and V. Vovchenko, J. Phys. G 42, 105102 (2015). 

[29] V. Vovchenko, D. V. Anchishkin, M. I. Gorenstein, and R. V. Poberezhnyuk, Phys. Rev. C 92 , 
054901 (2015). 

[30] M. Gazdzicki, M. I. Gorenstein, and M. Mackowiak-Pawlowska, Phys. Rev. C 88, 024907 (2013). 


16 


